Multi-variable special polynomials using an operator ordering method

Using an operator ordering method for some commutative superposition operators, we introduce two new multi-variable special polynomials and their generating functions, and present some new operator identities and integral formulas involving the two special polynomials. Instead of calculating complicated partial differential, we use the special polynomials and their generating functions to concisely address the normalization, photocount distributions and Wigner distributions of several quantum states that can be realized physically, the results of which provide real convenience for further investigating the properties and applications of these states.     

Sliding mode control for a diffuse interface tumor growth model coupling a Cahn–Hilliard equation with a reaction–diffusion equation

We consider the sliding mode control (SMC) problem for a diffuse interface tumor growth model coupling a Cahn–Hilliard equation with a reaction–diffusion equation perturbed by a maximal monotone nonlinearity. We prove existence and regularity of strong solutions and, under further assumptions, a uniqueness result. Then, we show that the chosen SMC law forces the system to reach within finite time a sliding manifold that we chose in order that the tumor phase remains constant in time.     

Numerical optimisation of chemotherapy dosage under antiangiogenic treatment in the presence of drug resistance

We consider a two-compartment model of chemotherapy resistant tumour growth under angiogenic signalling. Our model is based on the one proposed by Hahnfeldt et al. (1999), but we divide tumour cells into sensitive and resistant subpopulations. We study the influence of antiangiogenic treatment in combination with chemotherapy. The main goal is to investigate how sensitive are the theoretically optimal protocols to changes in parameters quantifying the interactions between tumour cells in the sensitive and resistant compartments, that is, the competition coefficients and mutation rates, and whether inclusion of an antiangiogenic treatment affects these results. Global existence and positivity of solutions and bifurcations (including bistability and hysteresis) with respect to the chemotherapy dose are studied. We assume that the antiangiogenic agents are supplied indefinitely and at a constant rate. Two optimisation problems are then considered. In the first problem a constant, indefinite chemotherapy dose is optimised to maximise the time needed for the tumour to reach a critical (fatal) volume. It is shown that maximum survival time is generally obtained for intermediate drug dose. Moreover, the competition coefficients have a more visible influence on survival time than the mutation rates. In the second problem, an optimal dosage over a short, 30-day time period, is found. A novel, explicit running penalty for drug resistance is included in the objective functional. It is concluded that, after an initial full-dose interval, an administration of intermediate dose is optimal over a broad range of parameters. Moreover, mutation rates play an important role in deciding which short-term protocol is optimal. These results are independent of whether antiangiogenic treatment is applied or not.     

A new non-perturbative approach in quantum mechanics for time-independent Schr?dinger equations

A new non-perturbative approach is proposed to solve time-independent Schr?dinger equations in quantum mechanics.It is based on the homotopy analysis method(HAM)that was developed by the author in 1992 for highly nonlinear equations and has been widely applied in many fields.Unlike perturbative methods,this HAM-based approach has nothing to do with small/large physical parameters.Besides,convergent series solution can be obtained even if the disturbance is far from the known status.A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be one thousand times larger than the possible superior limit of the perturbative approach.This HAM-based approach could provide us rigorous theoretical results in quantum mechanics,which can be directly compared with experimental data.Obviously,this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.     

Flow physics of normal and abnormal bioprosthetic aortic valves

Flow physics of transvalvular flows in the aorta with bioprosthetic valves are investigated using computational modelling. For the efficient simulations of flow-structure-interaction in transvalvular flows, a simplified, reduced degree of freedom valve model is employed with a sharp interface immersed boundary based incompressible flow solver. Simulations are performed for normal as well as abnormal valves with reduced leaflet motion that models the effect of early leaflet thrombosis. The structure of the aortic jet and the hemodynamic stresses on the aortic wall are analysed to understand the hemodynamic impacts and possible long-term clinical implications of sub-clinical, reduced leaflet motion. The simulation results have shown that the reduced leaflet motion tilts the direction of aortic jet and generates stronger flow separation and re-attachment on the aortic wall downstream from the reduced motion leaflets. The modified flow pattern increases the wall pressure fluctuation and average wall shear stress on the downstream aortic wall, and results in the asymmetric oscillatory shear index distributions, which may have long-term clinical implications such as aortic wall damage and remodelling.     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

砂岩方梁断裂过程区范围的实验确定方法

裂纹前端的断裂过程区是引起岩石非线性断裂及尺寸效应的主要原因。利用数字图像相关技术对砂岩开展了三点弯曲梁实验,获得观测区域高精度的全场位移和应变数据,根据断裂韧带区域水平位移和水平应变的分布特征,结合裂尖岩石颗粒变化的微观分析,提出采用裂纹尖端水平位移波动性和水平应变突变性所得到的波动系数和水平应变突变值,确定断裂过程区形状和临界尺寸的方法。结果表明:砂岩断裂过程区的形状为不规则的狭长带状区域,断裂过程区的临界长度为11~13mm,临界宽度为1.58~2.36mm。断裂过程区区域内形变在趋向裂尖时呈指数增加,但其单位区域内的形变增量呈波动状态。该方法能够更加准确判断岩石断裂过程区的范围,有助于分析岩石的非线性断裂特性。